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Abstract 



At HERA heavy quarks may contribute up to 30% of the structure function 
F2. The introduction of heavy quarks requires an extension of the DGLAP 
formalism. The effect of using different heavy flavour number schemes, and 
different approaches to the running of as, are compared using the ZEUS PDF 
fit formalism. The potential of including charm data in the fit is explored, using 
D* double differential cross-sections rather than the inclusive quantity Fg'^. 



Parton Density Function (PDF) determinations are usually global fits [1-4], which use inclusive 
cross-section data and structure function measurements from deep inelastic lepton hadron scattering 
(DIS) data as well as some other exclusive cross-ections. The kinematics of lepton hadron scattering is 
described in terms of the variables Q^, the invariant mass of the exchanged vector boson, Bjorken x, the 
fraction of the momentum of the incoming nucleon taken by the struck quark (in the quark-parton model), 
and y which measures the energy transfer between the lepton and hadron systems. The differential cross- 
section for the neutral current (NC) process is given in terms of the structure functions by 



where Y± = 1 it (1 — y)^. In the HERA kinematic range there is a sizeable contribution to the F2 
structure function from heavy quarks, particularly charm. Thus heavy quarks must be properly treated 
in the fomalism. Furthermore fitting data on charm production may help to give constraints on the gluon 
PDF at low-x. 

The most frequent approaches to the inclusion of heavy quarks within the conventional framework 
of QCD evolution using the DGLAP equations [5-8] areQ: 

• ZM-VFN (zero-mass variable flavour number schemes) in which the charm parton density c(x, Q"^) 
satisfies c{x,Q'^) = for < /i^ and n/ = 3 + 9{Q^ — ^1) in the splitting functions 
and (5 function. The threshold /i^, which is in the range < //^ < 4m^, is chosen so that 
F2 (x, Q"^) = 2e^xc(x, Q"^) gives a satisfactory description of the data. The advantage of this ap- 
proach is that the simplicity of the massless DGLAP equations is retained. The disadvantage is 
that the physical threshold = — 1) > is not treated correctly {W is the 7*5 CM 
energy). 

• FFN (fixed flavour number schemes) in which there is no charm parton density and all charmed 
quarks are generated by the BGF process. The advantage of the FENS scheme is that the threshold 
region is correctly handled, but the disadvantge is that large ln((5^/m^) terms appear and charm 
has to be treated ab initio in each hard process. 

• GM-VEN (general mass variable flavour number schemes), which aim to treat the threshold cor- 
rectly and absorb hi{Q'^ /ml) terms into a charm parton density at large Q^. There are differing 
versions of such schemes [9, 10] 

For the main ZEUS-S analysis [4], the heavy quark production scheme used was the general 
mass variable flavour number scheme of Roberts and Thorne (TR-VEN) [10, 11]. However we also 
investigated the use of the FFN for 3-flavours and the ZM-VFN. In Fig. [T] we compare the fit prediction 
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Fig. 1: ZEUS data on ^2'^ compared to predictions using tlie FFN (left) ZM-VFN(middle), TR-VFN (rigiit) sciiemes. In eacii 
c ase the fit parameters are kept the same (fi t ted using FFN) and only tlie scheme is ch a nged. 




Fig. 2: ZEUS data on Fl'^ compared to predictions using the FFN (left) ZM-VFN(middle), TR-VFN (right) schemes. In each 
case the fit parameters are refitted when the scheme is changed. 

for using each of these schemes to data on F^^^ from the ZEUS collaboration [12]. One can see the 
differences between the FFN and the ZM-VFN at threshold where the ZM-VFN is clearly inadequate. In 
this kinematic region the TR-VFN is more like the FFN. However, the TR-VFN scheme becomes more 
like the ZM-VFN scheme for >> m^. 

This comparison illustrates the effect of change in scheme when keeping the PDF parameters 
fixed. In practice one should refit the PDF parameters using the alternative schemes. The result of this is 
shown in Fig |2l The difference between the FFN and TR-VFN is not so marked. It is well known that 
these choices have some effect on the steepness of the gluon at very small-x, such that the zero-mass 
choice produces a slightly less steep gluon. In Fig [3] the differing shapes of the sea and the gluon PDFs 
for these different heavy quark schemes are illustrated. 

Figure |4] shows the ZEUS-S fit predictions for more recent Fg'^ data from ZEUS and HI [13, 14]. 



'charm production is described here b ut a similar formalism describes beauty pro d uction 




Fig. 3: The sea and gluon PDFs extracted from fits using the FFN (left) ZM-VFN(middle), TR-VFN (right) schemes. In each 
case the fit parameters are refitted when the scheme is changed 
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Fig. 4: Comparsion of ZEUS PDF fit predictions to recent charm data from ZEUS and HI on i^l'^ 
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Fig. 5: The gluon PDF and its fractional uncertainties at = lOGeV^, from a) the ZEUS-O PDF fit (left) and b) a smilar fit 
with Fi" data included (right). 

The scheme chosen was FFN for 3 flavours with the renormahsation and factorisation scale for light 
quarks both set to but the factorisation scale for heavy quarks set to + 4m^. The reason for 
these choices of scheme and scale is that these are the choices made in the programme HVQDIS [15-17] 
which was used to extract .Fl'^ from data on D* production. The scale choice does not make any signficant 
difference to the predictions (see later). 

Note that in Fig |4] the charm data are shown compared to the ZEUS-S PDF fit predictions but 
these data were not input to the fit. Including the ZEUS charm data [13] in the ZEUS-S PDF fit gives 
no visible improvement to PDF uncertainties. To investigate the potential of charm data to constrain the 
gluon PDF, we modified the ZEUS-S PDF fit as follows: all ZEUS inclusive neutral current and charged 
current cross-section data from HERA-I was included but no fixed target data; the parametrisation was 
modified to free the mid-x gluon parameter p;i{g) and the low-x valence parameter P2{u) = P2{d), 
however the d — u normalisation had to be fixed since there is no information on this without fixed 
taregt data. This fit is called the ZEUS-O fit. Fig. |5] compares the gluon PDF and its uncertainties as 
extracted from this ZEUS-O PDF fit with the those extracted from a similar- fit including the data. 
This illustrates that the charm data has the potential to constrain the gluon PDF uncertainties. Its lack of 
impact on the global fit may be because we are not using the charm data optimally. 

is a quantity extracted from D* cross-sections by quite a large extrapolation. It would be better 
to fit to those cross-sections directly. The evaluation of the theoretical predictions involves running the 
NLO programme HVQDIS for each iteration of the fit. However, one can shorten this process by using 
the same method as was used for the ZEUS -JETS fit [18]. The PDF independent subprocess cross- 



= 2 GeV' 
n HI 96-97 
■ HI prel. 
• ZEU.S 98-00 
o ZEUS 96-97 
ZEUSNLO 

:4r " 

|Tn!ir^-i 


' 4 GeV- 


' 7 GeV^ 


llGeV' 


18 GeV' 


1 1 1 

1 30 GeV' 






Iv - 


60 GeV' 


130 GeV^ 


500 GeV- 


•V 




s 




D* production 



Q? = 2 GeV^ 

' 1 1 1 < < 


' J% 

4 GeV^ 

' 1 1 1 < < 


- • ZEUS 98-00 

- ZEUS NLO Ag 
I — ZEUS NLO Ag 

+ Lund 

- A 
7 GeV^ • 




- f\ 

1 1 GeV^ 


: i\ 

18 GeV^ 


± 

- ri 

30 GeV^ 


60 GeV^ 


: Is 

130 GeV^ 


500 GeV^ ^ 



Fig. 6: Double differential cross-sections for D* production. The red lines show the predictions of the ZEUS-S-13 NLO 
PDF fit using the Petersen fragmentation function for the D* , whereas the blue lines show these predictions using the Lund 
fragmentation function. 



sections are output onto a grid, such that they can simply be multiplied by the PDFs at each iteration. 
The data used are the nine double differential cross-section measurements of (fa{D*)/dQ'^dy [13], see 
Fig. El 

There are further theoretical considerations to be accounted for when unputting D* cross-sections, 
as opposed to and inclusive quantity like F^^^, to a PDF fit. Since the grids are calculated using HVQDIS 
the fit must use the FFN scheme to be compatible. This means that we cannot use ZEUS high-Q^ 
data, since this scheme is not suitable at high-Q^. Hence we chose to use the ZEUS-S global fit, which 
incuded fixed target data, with a cut-off < SOOOGeV^. Furthermore, it has only recently become 
evident that since we are using the FFN scheme we must also treat the running of differently than 
in the VFN schemes. In these VFN schemes Ug is matched at flavour threholds [19], but the slope of 
Og is discontinuous at the flavour thresholds. For consistency with HVQDIS we must use a 3-flavour 
as which is continuous in Q^. This requires an equivalent value of as{Mz) = 0.105 in order to be 
consistent, at low Q^, with the results of using a value of as{Mz) = 0.118 in the usual VFN schemes. 
Such a 3-flavour as has also been used in specialised PDF fits of MRST (MRST2004F3) [20], which are 
used to make predictions for charm production. 

In Fig |7] we compare different heavy quark factorisation scales and different treatments of the 
running of ag for predictions of F^^. Fig.[8]makes the same comparision for We see that within 

the FFN scheme the choice of the heavy quark factorisation scale makes only a small difference at low 
Q^. The treatment of as gives larger differences. The FFN scheme and TR-VFN scheme differ for 
almost all if as runs as for the VFN schemes. However if a 3-flavour as is applied in the FFN 
schemes there is much better agreement of all schemes at higher Q^. 

We now return to consider inputting the D* cross-sections to the PDF fit. The ZEUS-S global 
fit formalism is used including all ZEUS inclusive neutral and charged current cross-section data from 
HERA-I and the fixed target data. The parametrisation was also modified to free the mid-x gluon pa- 
rameter P5{g) and the low-x valence parameter P2{u) = P2{d). This fit is called ZEUS-S-13. Figure|9] 
shows the difference in the gluon PDF uncertainties, before and afer the D* cross-sections were input to 
the ZEUS-S-13 global fit. Disappointingly the uncertainty on the gluon is NOT much improved. 

Should we have expected much improvement? There are two aspects of the fit which could be im- 
proved. The predictiond for the D* cross-sections have more uncertainties than just the PDF parametriza- 
tion. A further uncertainty is introduced in the choice of the c — > fragmentation The Petersen frag- 

^Note the predictions are always made by refitting PDF parameters for each scheme choice, not by simply changing the 
scheme with the same PDF parameters 
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Fig. 8: Comparison of predictions for F2 , from fits wliich use tlie TR-VFN sclieme and tlie FFN sclieme witli two different 
factorisation scales: on the left hand side the FFN schemes still use a VFN treatment of as, whereas on the right hand side a 
3-flavour Os is used. 




Fig. 9: The gluon PDF and its fractional uncertainties for various bins Left: before D* cross-section data are input to the 
ZEUS-S-13 fit. Right: after D* cross-section data are input to the ZEUS-S-13 fit 
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Fig. 10: Fractional uncertainties on Ddouble differential cross-sections for D* production, the red lines show the uncertainties 
onn these cross-sections deriving from the uncertainty on the gluon PDF in the ZEUS-S-13 fit, before including these D* data 
in the fit. 



mentation function was used for the fit predictions. However, looking back at Fig [6] we can see that the 
Lund fragmentation function seems to describe the data better. To best exploit the charm data in future 
we need to address such aspects of our model uncertainty. Secondly, Fig [TO] compares the fractional 
errors on the D* cross-sections with the uncertainty on the prediction for these quantities derived from 
the uncertainty on the gluon PDF in the ZEUS-S-13 PDF fit, before inputting the D* cross-sections. The 
data errors are larger than the present level of uncertainty. Thus we eagerly await the 5 -fold increase in 
statistics expected from HERA-II charm and beauty data. 
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